Chaotic Electron Motion in Superlattices. Quantum-Classical Correspondence 
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Abstract 

We investigate the classical-quantum correspondence for particle motion in a superlattice 
in the form of a 2D channel with periodic modulated boundaries. Its classical dynamics 
undergoes the generic transition to chaos of Hamiltonian systems as the amplitude of the 
modulation is increased. We show that for strong chaotic motion, the classical counterpart 
of the structure of eigenstates (SES) in energy space reveals an excellent agreement with the 
quantum one. This correspondence allows us to understand important features of the SES in 
terms of classical trajectories. We also show that for typical 2D modulated waveguides there 
exist, at any energy range, extremely localized eigenstates (in energy) which are practically 
unperturbed by the modulation. These states contribute to the strong fluctuations around 
the classical SES. The approach to the classical limit is discussed. 

PACS: 05.45+b, 03.20. 



A variety of important physical systems 
are periodic, either in space; such as crys- 
tals and superlattices; in time, hke periodi- 
cally forced mechanical devices; or in both, 
as in particle motion in a superlattice subject 
to an ac electric field. The analysis of such 
systems has lead to the understading of well 
known fundamental phenomena, e.g., energy 
band structure of solids and resonances. On 
the other hand, during the last decade much 
attention have been paid to to study the 
classical-quantum correspondence (CQC) of 
classically chaotic systems. In [1] we give 
references to some interesting physical time- 
periodic chaotic systems. One of the chal- 
lenges is to identify the ways in which classi- 
cal chaos may be manifested in the quantum 
regime [2]. There are well established tools 
to characterize a quantum classically chaotic 
system. The most commonly used relies on 
the analysis of their spectra (e.g., level spac- 
ing distribution) in the context of the ran- 
dom matrix theory conjecture [3]. See also 
[4] for other signatures of "quantum chaos". 

Here we explore the CQC of the particle 



motion in a superlattice in the form of a 2D 
electron waveguide with periodically modu- 
lated boundaries. Two possible experimen- 
tal realizations are: 1) a film whose thick- 
ness is a periodic function of one of the co- 
ordinates and 2) a periodically modulated 
mesoscopic electron channel in the ballistic 
regime (see [5,6] and references therein). We 
analyze the problem by studying the struc- 
ture of the eigenf unctions (SEF) and local 
density of states (LDOS) [7]. This novel ap- 
proach has been applied recently to a va- 
riety of systems [8]. The basic idea is as 
follows. Consider a non-integrable Hamilto- 
nian H = Ho+V, where V^is the non-integral 
perturbation to the integrable Hamiltonian 
Hq. Let the eigenstates of H and Ho be de- 
noted, respectively, by \E'" and cj)j. We can 
then form the matrix wf = |Cf (A;)p, where 
C° =< \E'° I (f)j >. The rows (columns) 
of show how a specific perturbed eigen- 
state I a > (unperturbed eigenstate >) 
is expanded in the unperturbed basis |/ > 
(perturbed basis a >). For our purposes 
it is essential, as will become clear below. 



to energy-order the unperturbed as well as 
the perturbed states; that is, E"'^^ > E"' 
and Ei_^^ > Ef, where E"' and E'f are the 
energy spectra of the perturbed and unper- 
turbed systems, respectively; aj = 1,2.... 
The structure of eigenstastes SES is defined 
as [7], 

W{Ef\E'^) = ^wf6{E^-E-'), (1) 

a' 

where ivf — wf /N with N as the num- 
ber of eigenstates \a' > in the vicinity of 
a given a. Eq. (1), seen as a function of 
the unperturbed energy Ef gives the struc- 
ture of eigenstates with total energy close 
to E'^. The local density of states LDOS, 
also known as the strength function, defined 
as t<;(E" I Ef) = E/'^?^(^" - E^'), is also 
used to study the quantum-classical corre- 
spondence [6,8,9]. Here, for lack of space we 
shall not discuss it. 

Now, the reason these two quantities, SES 
and LDOS, help us tackle the quantum- 
classical correspondence is that they have 
well defined classical counterparts [7] . Since 
=< I 0; > as a function of / is the 
projection of the perturbed state onto the 
states of the unperturbed system, the clas- 
sical counterpart of =| p as a func- 
tion of energy Ei can be defined as the pro- 
jection of the total Hamiltonian H onto the 
unperturbed one H^, where H = + V 
with V the perturbation. This can be nu- 
merically done by substituting the trajecto- 
ries = {x(t),y(t),p^(t),py(t)) generated 
by H with energy E into Since the 

unperturbed energy E^{t) along these tra- 
jectories varies in time, it fills the so-called 
energy shell characterized by its width AE. 
For chaotic total Hamiltonians H, the classi- 
cal analog Wci{E'^ \ E) of the quantum SES 
is the probability distribution constructed 
from E^{t). See below for details. 

Similarly, the classical analog of LDOS can 
be obtained by performing the inverse opera- 
tion, namely, projecting the dynamics of 
onto the Hamiltonian function H. 



This approach assumes that the Hamilto- 
nian could be separated into unperturbed 
and perturbed parts. For billiards, like our 
system, the non-integrability comes from the 
boundary conditions; the Hamiltonian op- 
erator, being simply the kinetic energy, is 
the same for both perturbed (periodically 
modulated) and unperturbed (flat) channel. 
To overcome this problem, a transformation 
to curvilinear coordinates can be performed 
such that both, perturbed and unperturbed 
systems have the same boundary conditions. 
Consequently the effects of the boundary ap- 
pear in the new coordinates as an interaction 
potential between the two degrees of freedom 
(see details in [9]). 

Since the Hamiltonian is periodic, the en- 
ergy eigenstates satisfy Bloch's theorem. Al- 
though most of the results will be indepen- 
dent of the details of the profile, we need to 
specify a particular one. As in [5,6,9] we de- 
fine the top boundary hj y = d + a cos(x) 
and the bottom boundary hy y = 0. The 
first classical dynamical studies of this sys- 
tem appear in [10,11]. The finite length ver- 
sion this system was analyzed in [12] as a 
model of a mesoscopic electron waveguide, 
where a transport signature of chaos in the 
behavior of resistivity was established. Ket- 
zmerick [13] also used it to study certain 
quantum transport properties through bal- 
hstic cavities. Moreover, the analysis of the 
band-energy spectrum for an infinitely long 
rippled channel [5] in the cases of mixed and 
global classical chaos, yields insight into the 
universal features of electronic band struc- 
tures of real crystals [14]. See also [15]. Re- 
sults obtained for this particular system are 
applicable to a large class of systems, namely 
non-degenerate, non-integrable Hamiltoni- 
ans. We hmit this report to the case of global 
classical chaos which occurs in wide {d > 27r) 
channels [5]. The geometrical parameters 
used here arc {a/2n,d/2n) = (0.12,0.5). 
Fig. 1 shows the unperturbed energy E'^{t) 
obtained as explained above. The distribu- 
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tion of energies on the right of Fig. 1 is, by 
definition, the classical counterpart of SES, 
obtained by averaging E^{t) over time. 
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FIG. 1 a) E^{t) as a function of time (in arbitrary 
units) for £ = 1. b) Classical SES Wci{E^ \ E). 




FIG. 2 Structure of eigenstates. Quantum SES 
(black line); Classical SES (white line). The 
average for the SES is over 1900 <a< 1990. Inset: 
classical SES for (a/27r, rf/27r) = (0.06, 1) (dotted) 
and (a/27r,d/27r) = (0.12,0.5) (continuous). 

Fig. 2 compares the quantum and classical 
SES. We see that the classical SES agrees 
well with the average shape of quantum SES. 
I.e., the classical SES predicts: 1) the ap- 
pearance of three prominent peaks at spec- 
ified values, 2) the width of the distribu- 
tion (the energy shell); and 3) its assym- 
metry about 1. The inset compares classi- 
cal SES for two channels defined by different 
geometrical parameters but both displaying 
global chaos in phase space. Thus, the clas- 
sical and quantum SES can detect dynam- 
ical differences not revealed by the standard 
tools, e.g., Poincare maps, and thus serves 



to complement the characterization of dy- 
namical systems. In fact, detailed analysis 
[9] demonstrates that the left (right) side 
peak is formed by trajectories dwelling near 
the unstable (stable) fixed period-one fixed 
point, whereas the central peak is due to 
grazing trajectories. 

The strong fluctuations of the quantum 
SES are due to: 1) the fact that the system 
is still not deep in the semiclassical regime 
(although the level numbers are high) and 
2) the existence of extremely localized and 
sparse (in energy space) states. See below 
and [9] for details. 
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FIG. 3 Typical eigenfunctions. a) a=1984, b) 
a=1985, and c) a=1986. 

Fig. 3 shows 3 consecutive eigenstates in 
the energy-ordered unperturbed basis, illus- 
trating the kind of eigenstates that can occur 
typically: extended, sparse, and localized in 
the energy representation. A quantitative 
characterization is given by various localiza- 
tion measures, such as the entropy localiza- 
tion length Ih, the inverse participation ratio 
lipri and the mean square root l„. See [6] for 
their definitions. In Fig. 4 we present In as 
a function of the perturbed state. It shows 
wild fluctuations, due to the existence of the 
localized and sparse states mentioned above. 
The extremely localized states are those that 
are practically the same as the unperturbed 
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(flat channel) states with transversal mode 
number m = 1. It can be shown that the 
number of these states A^rre=i ; relative to the 
total number of states N{E), up to energy 
E, is given by Nm=i/N(E) = A/ttVE. 
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FIG. 4 (a) Entropy localization length Ih as a 
function of a. (b) Nm=i{E)/N{E) as a function of 
E (dots) for the first 2000 states. The solid line is 
the analytical estimate A/tt^/E. 

The fact that Nm=i/N{E) decays as 1/y/E 
implies that in the strict classical limit 
the fluctuations, observed in Fig. 2 for 
the quantum SES, will vanish, however, 
the approach to this limit is surprisingly 
slow. Furthermore, it can be shown that 
Njn=i = 27r\/E/d, where d is the width of 
the channel, in units of the period of the 
ripple. This indicates that localized states 
continue to appear at all energies, even in 
the classical limit!. That is, one can always 
flnd states that remain practically unper- 
turbed by the modulation. The importance 
of these extremely localized states in finite 
electron waveguides with modulated or cor- 
rugated boundaries is addressed in [16]. 
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